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METHOD AND APPARATUS PROVIDING IMPROVED 
POSITION ESTIMATE BASED ON AN INITIAL COARSE 

POSITION ESTIMATE 

BACKGROUND 

Field 

[1001] The present invention relates generally to location of a device, and more 
particularly to techniques for providing an improved (i.e., more accurate) estimate of the 
location of the device based on an initial coarse position estimate. 

Related Applications 

[1002] This application claims priority to a U.S. Provisional Application Serial No. 
60/273,570, filed on March 5, 2001, and assigned to the assignee of the present 
application. 

Background 

[1003] A common means by which to locate a device is to determine the amount of 
time required for signals transmitted from multiple sources at known locations to reach 
a receiver within the device to be located. One system that provides signals from a 
plurality of transmitters of known locations is the well-known Global Positioning 
Satellite (GPS) system. Satellites in the GPS system are placed in precise orbits 
according to a GPS master plan. The position of the GPS satellites can be identified by 
a number of different sets of information^ some more accurate than others. 
[1004] GPS satellites transmit a set of information, referred to as "Almanac'% whidi 
includes less accurate information regarding the location of the satellites in the 
"constellation". Ground stations continuously monitor the GPS satellites to observe 
their variations in orbit. Once the satellite positions have been measured, the 
information is relayed back to the satellites. The satellites then transmit another set of 
information, refeired to as '"Ephemeris", which includes a higher accuracy version of 
the satellite orbits. Each satellite transmits the Almanac information for all satellites but 
the Ephemeris information only for itself. 
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[1005] A GPS receiver can receive and/or store an almanac that indicates where 
each of a number of satellites is located in the sl^ at a given time. A more accurate 
determination of the location of a GPS satellite can be made based on the Ephemeris 
and the time of day at which this information is available. 

[1006] Both the Almanac and Ephemeris are valid for a limited amount of time. 
Hie Almanac information is considered to be accurate to ^ptx)ximately 3 kilometers fdr 
approximately one week from the time the Almanac is transmitted. The Ephemeris 
provides infomiation regarding the satellite orbit with an accuracy of approximately 1 
meter for approximately 2 hours. The error in both the Almanac and Ephemeris grows 
as the information ages. Accordingly, the location of the satellites based on this 
information is less and less accurate as the Almanac and Ephemeris age, unless updated 
information is received at appropriate intervals in time. 

[1007] Without accurate information regarding the location of the satellites, the 
estimated location of a device, which may have been determined based on the receipt of 
signals transmitted from the satellites, may be inaccurate. Accurate information may be 
attained by receiving updates (continually or as necessary) from the satellites or from an 
alternative source. The alternative source may be a base station or position determining 
equipment (PDE) in a wireless communication system, either of which may have a GPS 
receiver capable of receiving the required information from the GPS satellites. 
However, valuable resources would be consumed for the device to be located to attain 
die information at regular intervals. In particular, power is required to transmit and 
receive the information, and bandwidth is required to transmit the information from a 
remote source to the device. 

[1008] There is therefore a need in the art for techniques to accurately determine the 
position of a device with minimal expenditure of resources. This need is particularly 
acute when position determination is based on transmitters (e.g., GPS satellites) that 
move over time and when accurate locations of such transmitters are known only upon 
receiving updates from the transmitters or a source remote from the satellites and the 
device whose position is to be determined. 

SUMMARY 

[1009] Techniques to accurately determine the location of a receiver device based 
on an initial coarse position estimate, which may have been derived using less accurate 
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infonnation regarding the location of the transmitters (e.g., an old Almanac or old 
Ephemeris for the GPS satellites) are disclosed herein. In one embodiment of the 
disclosed method and apparatus, corrections to the coarse position estimate of a receiver 
device are made based on knowledge of: (1) information providing a relatively more 
accurate location of the transmitters; and (2) information providing a relatively less 
accurate location of the transmitters (e.g., the old Almanac and/or Ephemeris) used to 
derive the coarse position estimate. The corrections may be performed based on various 
correction algorithms, one of which is an iterative algorithm described in further detail 
below. 

[1010] In accordance with one embodiment of the disclosed method and apparatus, 
a coarse position estimate for the device is received. The coarse position estimate may 
have been derived based on initial O&ss accurate) estimates of the position of a plurality 
of transmitters, such as GPS satellites. Revised (more accurate) position estimates for 
the transmitters are also received. The revised position estimate for the device is then 
initialized (e.g., to the coarse position estimate). An update vector is next computed 
based on the initial and revised position estimates for the device and the initial and 
revised position estimates for the transmitters. The revised position estimate for the 
device is then updated based on the update vector. The computation for the update 
vector and the updating of the revised position estimate for the device can be repeated a 
number of times (e.g., until the magnitude of the update vector is within a particular 
threshold) to achieve a more and more accurate estimate of the actual position of the 
device. 

[1011] In accordance with one embodiment of the disclosed method and apparatus, 
the update vector is computed by performing the steps of: (1) calculating a set of 
measurements based on the revised position estimate for the device and the revised 
position estimates for the transmitters; (2) deriving an int^mediate position estimate for 
the device based on the set of measurements and the initial position estimates for the 
transmitters; (3) determining a first geometry matrix based on the initial position 
estimates for the transmitters and the intermediate position estimate for the device; (4) 
determining a second geometry matrix based on the revised position estimates for the 
transmitters and the revised position estimate for the device; and (5) computing the 
update vector based on the intermediate and revised position estimates for the device 
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and the first and second geometry matrices. These steps are described in further detail 
below. 

[1012] The disclosed method and apparatus further provides otiier aspects, 
embodiments, and features, as described in fuither detail below. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[1013] The features, nature, and advantages of the present invention will become 
more apparent from the detailed description set forth below when taken in conjunction 
with the drawings in which like reference characters identify like elements. 
[1014] FIG. lA is a simplified illustration of a system capable of inq>lementing 
various aspects of the invention. 

[1015] FIG. IB is a diagram that illustrates the error in the location of the 

transmitters based on information available to the receiver device. 

[1016] FIG. 2 is a simplified block diagram of the receiver device and the position 

identifying device in accordance an embodiment of the invention. 

[1017] FIG. 3 is a flow diagram illustrating a process performed by the position 

identifying device to derive a more accurate position estimate for the receive device 

based on a coarse position estimate previously derived by the receiver device. 

[1018] FIG. 4 is a flow diagram of the processing performed for the iterative 

algorithm, in accordance with an embodiment of the invention. 

DETAILED DESCRIPTION 
[1019] FIG. lA is a simplified illustration of a system capable of implementing 
various aspects of the invention. A receiver device 110, the position of which is to be 
determined, receives signals transmitted from a plurality of the transmitters 130a-130d 
(referenced collectively using reference number "130") via an antenna. The receiver 
device 110 further communicates with a position identifying device 120, which assists 
in determining the position of the receiver device, as described in further detail below. 
[1020] In one embodiment, the receiver device 110 is a cellular telephone capable 
of receiving signals from the plurality of the transmitters 130. However, the receiver 
device 110 may be any device capable of determining the arrival times of received 
signals with respect to a reference time. For example, the receiver device 110 may be a 
computer terminal having a wireless modem, a stand-alone GPS receiver, a receiver 
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capable of receiving signals firom ground-based transmitters, or any other receiver. The 
transmitters 130 may be any type of transmitter having locations that are known or 
which can be ascertained. In an embodiment, the transmitters 130 are satellites of a 
Global Positioning Satellite (GPS) system. In alternative embodiments of the disclosed 
method and apparatus, the transmitters 130 may be ground-based transmitters (e.g., base 
stations of a wireless communication system), or a combination of GPS and ground- 
based transmitters. 

[1021] The receiver device 110 estimates its position based on the received signals 
and information indicative of the location of the transmitters 130 from which the 
received signals were transmitted. Due to various factors, the receiver device 110 may 
not have current or accurate information regarding the actual location of the transmitters 
130. In such case, the position estimate made by, the receiver device 110 of its location 
is coarse and may not have the desired degree of accuracy. The initial coarse position 
estimate made using less accurate transmitter locations may thereafter be "corrected" to 
provide a more accurate position estimate of the receiver device using the techniques 
described herein. 

[1022] The transmitters 130 can be used as reference points to determine the 
location of the receiver device 110. By accurately measuring the distances from three 
transmitters 130, the receiver device 110 can "triangulate" its position. The receiver 
device 110 determines distance by measuring the time required for a signal to travel 
from a transmitter 130 to the receiver device 1 10. By knowing the time the signal is 
transmitted from the transmitter 130 and observing the time the signal is received by the 
receiver device 1 10 (based on its internal clock), the travel time of the signal can be 
determined. However, the exact amount of time between transmission and reception 
typically cannot be determined, for example, due to offsets in the two clocks at the 
transmitter 130 and the receiver device 110. Thus, a "pseudo-range'* is typically 
computed based on the difference between a "reference" time and the time Aat die 
signal is received. The reference time may be any time, as long as the reference time is 
common to all pseudo-range measurements being made (or the pseudo-range 
measurements can be adjusted to compensate for differences in the reference times 
used). 

[1023] In FIG. lA, the transmitters 130a through 130d are shown at positions 
estimated by the receiver device 110 based on less accurate (e.g., not up-to-date) 
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infonnation. These positions may be different from the actual positions of the 
transmitter 130. In the instance in which the transmitters 130 are satellites, such as GPS 
satellites, the position of the satellites can be identified by the Almanac and Eph^neris. 
[1024] An accurate determination of the location of a GPS satellite can be made if 
the current Ephemeris information and the time of day are available. However, the 
Ephemeris information is valid for a limited amount of time (e.g., two hours). A less 
accurate determination of the GPS satellite location can be made if current Almanac 
information and the time of day are available. However, the Almanac is also valid for a 
limited amount of time (e.g., one week). 

[1025] Aspects of the invention provide techniques to accurately determine the 
location of a receiver device based on an initial coarse position estimate, which may 
have been derived using less accurate information regarding the location of the 
transmitters (e.g., an old Ahnanac or old Ephemeris for the GPS satellites). It should be 
understood that the initial estimate may be made based upon completely outdated 
Almanac or Ephemeris information. 

[1026] Since the pseudo-ranges and locations of the transmitters are used to 
determine the location of the receiver device, errors or inaccuracies in the locations of 
the transmitters translate into errors or inaccuracies in the position estimate for the 
receiver device. Hence, a "coarse" position estimate is derived for the receiver device 
based on pseudo-ranges to transmitters 130, the position of which may be known with 
relatively poor accuracy. 

[1027] In accordance with one embodiment of the disclosed method and apparatus, 
corrections to the coarse position estimate of receiver device are made based on: (1) 
more accurate information regarding the location of the transmitters 130; and (2) 
information regarding the amount of error in the estimates of the location of the 
transmitters 130 used to derive the coarse position estimate. Techniques to perform the 
corrections are described in further detail below. In one embodiment, the corrections 
are made at the position identifying device 120, which may be remotely located with 
respect to the receiver device 110. However, in other embodiments of the disclosed 
method and apparatus, the position identifying device 120 may be co-located with the 
receiver device 110. 

[1028] FIG. 2 is a simplified block diagram of the receiver device 110 and the 
position identifying device 120 in accordance with one embodiment of the disclosed 
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method and ^paratus. The receiver device 110 may be a component of a remote 
terminal in a wireless conmiunication system, such as a cellular telephone or a computer 
with a wireless modem. Alternatively, the receiver device 110 may be a stand-alone 
position determining unit, such as a stand-alone GPS receiver. 

[1029] The position identifying device 120 may be a component of a base 
transceiver subsystem (BTS), a base station controll^ (BSC), or a mobile switx^ng 
controller (MSG) in a wireless communication system, or may be a device that couples 
to one of these elements of the wireless communication system. For example, the 
position identifying device 120 may be a Position Determining Equipment (PDE) in a 
wireless communication system. Alternatively, the position identifying device 120 may 
be a component of, and co-located with, a remote terminal or some other device that 
includes the receiver device 110. Alternatively, the position identifying device 120 may 
be a stand alone component. 

[1030] For simplicity, the receiver device 110 and the position identifying device 
120 are shown in FIG, 2 as being directly coupled. However, one or more additidhal 
elements (e.g., a BTS, a BSC, and so on) may be coupled between the receiver device 
110 and the position identifying device 120. The position identifying device 120 and 
the receiver device 110 may also be co-located. This may be the case in a system in 
which accurate information about the location of the transmitters 130 is not available 
when the measurements are made, but is available at some time in the future. A coarse 
position estimate for the receiver device 110 may be derived immediately upon making 
the pseudo-range (or actual range) measurements. The coarse position estimate may be 
stored until information regarding the more accurate location of the transmitters 130 
becomes available, after which a more accurate position estimate for the receiver device 
may be derived. 

[1031] In the embodiment shown in HG. 2, the receiver device 110 includes a 
receivo- 210, a processing circuit 212, a memory 214, and a communications port 216, 
and further douples to an antenna 208.^ The antenna 208 receives signals Ax>m the 
transmitters 130 (shown in FIG. lA) and couples the received signals to the receiver 
210. The receiver 210 includes circuitry (e.g., radio finequency (RF) processing circuitry 
and/or other receiving circuitry) necessary to process the received signals to derive 
information used for position determination. Such information may include timing 
information and so on. The information is provided to the processing circuit 212, which 
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perfonns position detmnining calculations. Such calculations may include calculations 
to derive pseudo-ranges to the tiransmitters 130. Alternatively or additionally, the 
calculations may include those required to derive an initial coarse estimate of the 
position of the receiver device 110. The coarse position estimate may have a relatively 
large error due to large errors in the position estimates of the transmitters 130. 
[1032] The receiver device 110 provides the results of the calculations and other 
pertinent information via a communication port 216 to the position identifying device 
120. The p^nent information may include the time the received signals were 
measured, the information used to perform the calculations (e.g., the particular Almanac 
used to estimate the position of the transmitters 130), and possibly other information. 
The communication port 216 provides an interface to the position identifying device 
120. For a receiver device disposed within a remote terminal of a wireless 
communication system, the communication port 216 may support communication via a 
wireless link. 

[1033] In the embodiment shown in FIG. 2, the. position identifying device 120 
includes a conmaunication port 220, a processing circuit 222, and a memory 224. The 
communication port 220 receives the calculation results and pertinent information from 
the receiver device 110 and forwards the information to the processing circuit 222. The 
processing circuit 222 receives the calculation results and pertinent information, and 
further receives a copy of the "less accurate" information that was used by the receiver 
device 110 to perform the calculations. The less accurate information may be the 
Almanac (and the time at which the Almanac was valid), the Ephemeris (and the time at 
which the Ephemeris was valid), and/or any other information that the receiver device 
110 may have used to estimate the location of the transmitters 130. 
[1034] Tlie less accurate information may be provided by the receiver device 110 
together with the calculation results, e.g., as part of the pertinent information. In certain 
embodiments, the less accurate information may not be explicitly provided, and other 
indicators may be used to deduce the information that was used For example, the time 
for which the Almanac or Ephemeris is valid, the time at which the received signals 
were measured at the receiver device 1 10, the time at which the coarse position estimate 
was sent to the position identifying device 120, or some other information (which may 
be sent as part of the pertinent information), may be sufficient to allow the position 
identifying device 120 to correctly deduce the transmitter position estimates that were 
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used to derive the coarse position estimate. In some other embodiments, the position 
identifying device 120 may be responsible for sending to the receiver device 110 the 
information used to derive the coarse position estimate. For these embodiments, diere 
would be no need for the receiver device 110 to inform the position identifying device 
120 what information was used. 

[1035] In certain embodiments, the processing circuit 212 may communicate 
directly with the position identifying device 120 without the need for a separate 
communication port Tliis may be the case in instances in which the position 
identif3dng device 120 and the receiver device 110 are co-located, or may even be 
possible in certain other instances in which the position identifying device and receiver 
device are not co-located. 

[1036] In one embodiment, the position identifying device 120 receives die less 
accurate information prior to receipt of the calculation results from the receiver device 
110. In some embodiments, the less accurate information used by the receiver device 
110 may be provided to the position identifying device 120 by a source other than the 
receiver device 110 (e.g., over a communication link that is not shown in FIG. 2 for 
simplicity). For example, the less accurate information may be received directly by the 
position identifying device 120 from the transmitters 130. Alternatively, the less 
accurate information may be received by the position identifying device 120 from a 
source that is distinct from any of the components shown in FIG. lA, such as a 
component of a base station in a wireless communication system. 

[1037] In one embodiment, the position identifying device 120 may have multiple 
sets of information (e.g., several versions of the Almanac), any one of which may have 
been used by the receiver device 110 to perform the position determining calculations. 
In this case, the receiver device 110 may need to provide additional information to the 
position identifying device 120 to indicate what information, from among those 
available to the position identifying device 120, was specifically used by the receiver 
device 1 10 to perform the position determining calculation. 

[1038] FIG. 3 is a flow diagram illustrating a process performed by the position 
identifying device 120 to derive a more accurate position estimate for the receiver 
device 110 based on a coarse position estimate previously derived by the receiver 
device. The position identifying device 120 initially receives from the receiver device 
110 the results of the position determination calculations (e.g., the coarse position 
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estimate u"), at step 312. The position identifying device 120 also receives 
information indicative of the initial position estimates of the transmitters 130 (i.e., 
the less accurate transmitter position estimates), which were used to derive the coarse 
position estimate, at step 314. This information may indicate the version of the 
Almanac or Ephemeris that was used to derive the coarse position estimate. The initial 
position estimates for the transmitters 130, as determined by the receiver device 1 10, are 
less accurate and correspond to the location where the receiver device 110 assmned the 
transmitters to be located when making the position determining calculations, in an 
embodiment, the position identifying device 120 further receives the time at which the 
pseudo-range measurements were taken, also at step 314. 

[1039] The position identifying device 120 then detemfiines more accurate estimates 
sf of the location of the transmitters 130, at step 316. These more accurate transmitter 
position estimates may be made based on an Almanac and/or Ephemeris that is more 
accurate for the time at which the pseudo-range measurements were made by the 
receiver device 110. In one embodiment, the position identifying device 120 maintains 
a log of the Almanac and Ephemeris transmitted by the satellites of the GPS 
constellation. Such a log allows the position identifying device 120 to use the most 
accurate Almanac and Ephemeris to correct the coarse position estimate received firom 
the receiver device 110, as described in further detail below. 

[1040] Once the position identifying device 120 has information regarding the 
initial less accurate position estimates and the more accurate position estimates for the 
transmitters 130, the position identifying device 120 then corrects the coarse position 
estimate for the receiver device 110 in accordance with an algorithm described below, at 
step 318. 

[1041] FIG. IB is a diagram that illustrates the error in the location of the 
transmitters 130 based on information available to the receiver device 110. FIG. IB 
shows transmitters 132a through 132d at locations that represent the more accurate 
estimates of the actual location of the transmitters. FIG. IB also shows the transmitters 
130a through 130d at locations that represent less accurate estimates of the location of 
the transmitters. Transmitter 132 (shown with dashed lines) and the transmitter 130 
(shown with solid lines) correspond to two position estimates for the same transmitter. 
However, the location of each transmitter is perceived to be different depending upon 
whether the location of the transmitter has been determined using the more or less 
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accurate infonnation (i.e., using Almanac and Ephemeris that is, or is not, accurate for 
the time at which the pseudo-range measurement was taken). The locations of the 
transmitters 130 were used by the receiver device 110 to perforai the position 
determining calculations to derive the coarse position estimate for the receiver device 
110. 

[1042] The position identifying device 120 performs corrections on the coarse 
position estimate from the receiver device 110 to provide a more accurate position 
estimate for the receiver device. The corrections on the coarse position estimate may be 
perfonned based on various algorithms including a "linearized" algorithm, an "iterative 
algorithm", and possibly others. The linearized algorithm is described in detail in U.S 
Patent Application Serial No. 09/773,207, entifled "METHOD AND APPARATUS 
FOR DETERMINING LOCATION USING A COARSE POSIHON ESTIMATE," 
filed January 30, 2001, assigned to the assignee of the present application. The iterative 
algorithm is described below. 

[1043] The iterative algorithm can be used to perform corrections on a coarse 
position estimate for the receiver device 110 based on: (1) knowledge of the less 
accurate estimates of the location of the transmitters 130 (e.g., GPS satellites) used to 
derive the coarse position estimate for the receiver device 110, and (2) a knowledge of 
the more accurate estimates of the actual location of the transmitters 130 at the time the 
pseudo-range measurements were made. The more accurate transmitter position 
estimates may be determined based on knowledge of the specific time at which the 
pseudo-range measurements were taken. 

[1044] It should be understood that the coarse position estimate is made at either a 
location or a time when the more accurate estimates of the location of the transmitters 
130 are not available. Corrections wiU typically be performed at either a later time or 
another location or both, depending upon when and where the more accurate estimates 
of the position of the transmitters 130 become available. 

[1045] In one embodiment of the disclosed method and apparatus, it is assumed that 
there is a one-to-one mapping between a position location solution (i.e., the calculated 
position for the receiver device 110) and a set of transmitter position estimates. This 
assumption is easily satisfied in the case where there are only four pseudo-range 
measurements. 
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[1046] For clarity, the iterative algorithm is described specifically for a case in 
which the transmitters 130 are GPS satellites. However, the iterative algorithm may 
also be used with any types of transmitter in which inaccurate transmitter location 
information is available at a particular time or place, and more accurate transmitter 
location information is available at some later time and/or other location. 
[1047] The location s° for the.i-th satellite derived based on the Almanac 
information and the clock bias may be expressed as: 



Clock bias is defined as the difference between the time indicated by the local clock in 
the receiver device 110 and the time indicated by the clock in the transmitter 130. The 
location, will typically have a relatively large error due to inaccuracy of the 

Almanac used by the receiver device 1 10 to derive the location, s/* . 

[1048] The location, J* calculated for the i-th satellite based on the Ephemeris 

information and the clock bias may be expressed as: 



The location, s' will typically have a relatively small error due to the use of the more 

accurate Ephemeris information. The coarse position estimate, u" derived by the 
receiver device 110 using the relatively less accurate Almanac information may be 
expressed as: 



A more accurate position estimate u' of the receiver device 110 derived by the receiver 
device 110 using the relatively more accurate transmitter location estimates may be 
expressed as: 




r = ^\ y\ z'l and b',. . 
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Initially, more accurate infonnation may not be available. Therefore, the value of the 
more accurate position estimate may be set to . The position estimate is 
thereafter updated as more infonnation is obtained. 

[1049] A Range estimate associated with each satellite may be obtained based on 
the following expression: 



These range estimates are essentially pseudo-range measurements that would have been 
determined by a receiving device located at receiving signals &om transmitters 
located at the locations indicated by the more accurate infoimation. It can be seen that 
the measurement, nti of equation (1) is the distance between the more accurate position 
estimate u' of the receiver device 110 and the more accurate position estimate of 

one of the transmitters 130, taking into account the clock bias. Based on tiie range 
values of m,* calculated to each satellite, a least mean square (IMS) computation can be 
performed to determine the position = [x,' y^ ] that the receiver device 1 10 

would have calculated if it had been at the location and had assumed the transmitters 
130 to be at the position indicated by the less accurate position estimates, J/* for the 

transmitters 130 (i.e., the positions indicated by the Almanac information). Based on 
this position estimate it and the knowledge of the partial derivative equation between 

and M " , an update vector (SI can be determined. 
[1050] A geometry matrix for the less accurate transmitter position estimates 

^1" (e.g.. Almanac-derived satellite locations) and the corresponding location it 
derived from the previous LMS computation can be defined as: 



Eq(l) 




= K-^i 



K-^si yu-ys2 Zg-zl 



1 



Eq(2) 



M 



M 



M 



M 
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Similarly, a geometry matrix for the more accurate transmitter position estimates 
s' (e.g., the Ephemeris-derived satellite locations) and the more accurate receiver 
device position estimate can be defined as: 



X,, 



M 



■ys2 



e e 



M 



M 



jiT-^l ^-^1 ^-s:\ 



M 



-1 



Eq(3) 



[1051] A final geometry matrix H, can then be expressed as: 
H,={hI -R'^ H^)^ Hi -H^ , 



Eq(4) 



where R is the covariance matrix of the measurements, i.e., R = E[»r - J where e\ ] 
represents calculating the statistical expectation of the quantity in the brackets and m is 
the vector of measurements determined in equation (1). The update vector du can then 
be expressed as: 



Eq(5) 



[1052] The more accurate estimate of the receiver device's position obtained 
using the more accurate transmitter position estimates sf (e.g., the Ephemeris-derived 
satellite position), can be expressed as: 



Eq(6) 



Equations (1) through (6) can be applied iteratively a number of times to derive a more 
and more accurate position estimate for the receiver device 1 10. 

[1053] FIG. 4 is a flow diagram of the processing performed for the iterative 
algorithm, in accordance with an embodiment of the invention. The coarse position 
estimate derived by the receiver device 110 based on less accurate position 
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estimates J" of the transmitters 130 is iirst received (as described in FIG. 3). The more 
accurate position estimate for the receiver device 110 is initialized to the coarse 
position estimate m " , at step 412. Pseudo-range measurements m,- are then determined, 
at step 414. These measurements are the ones that would have been obtained if the 
position estimate were the solution. The measurements m/ can be derived based on 
equation (1) and can be detemained for all transmitters previously used to derive the 
coarse position estimate , 

[1054] At step 416, an LMS computation is performed based on the measurements 
nii determined in step 414 and the less accurate transmitter position estimates s° . The 
LMS computation can be achieved in a manner known in the art. The resulting solution 
from the LMS computation is denoted as u\ which is the solution that the receiver 
device 1 10 would have obtained if it had been at location u' . 

[1055] The geometry matrix H, is then computed in accordance with equations (2), 
(3) and (4). at step 418. The update vector dSF can then be determined based on the 
geometry matrix and the difference in position estimates (i.e., ir"-«') and in 

accordance with equation (5), at step 420. The more accurate position estimate u' is 
then updated with the update vector du , at step 422. 

[1056] At step 424, a determination is made whether or not the magnitude of the 
update vector du is less than a particular threshold (e.g., \du\ < 1). If the magnitude of 

the update vector du is not less than the threshold, the process returns to step 414 and 
the measurements m,- of the transmitters 130 are again determined based on the recently 
updated more accurate position estimate for the receiver device 110. Steps 414 
through 424 are repeated as many times as necessary until the magnitude of the update 
vector du is less than the threshold. Once the threshold is satisfied, the position 
estimate is returned as the more accurate location estimate for the receiver device 
1 10. The processing then terminates. 

[1057] The techniques described herein allow the location of a receiver device to be 
accurately determined based on a coarse location estimate for the receiver device, which 
may have been derived based on less accurate position estimates for the transmitters. In 
this manner, the information descriptive of the location of the transmitters (e.g., the 
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Almanac) may be loaded infrequently into the receiver device, which may allovf for 
conservation of resources. 

[1058] Hie techniques described herein may also be used advantageously for 
detennining the position of a receiver device operating in an analog mode, which is 
characterized by the time stamp of the measurements not being known with sufficient 
accuracy (e.g., more than 5 msec of time stamp error). 

[1059] In accordance with an aspect of the invention, measurements for five or 
more transmitters (e.g., GPS satellites) may be used to solve for x, y, z, and clock bias 
as well as the time-stamp error. The correction algorithm described above can be made 
applicable for cases where the time-stamp reported by the receiver device is not 
accurate. 

[1060] The iterative algorithm described above can be implemented for position 
determining calculations based on five or more measurements. The computations 
described above in equations (1) through (6) and the flow diagram shown in FIG. 4 can 
be applied in the manner described above. However, instead of the LMS algorithm 
normally used for four measurements, a modified LMS algorithm is implemented to 
operate on five unknowns to derive the solution (i.e., the position estimate for the 
receiver device). Instead of a dimensionality of four for geometry matrices Ha and He 
respectively shown in equations (2) and (3), these matrices would have a dimensionality 
corresponding to the number of measurements used. 

[1061] When five measurements are used for the position determining calculations, 
the geometry matrices Ha,5 can be expressed, as: 





' a 

y,, - ysi 


' a 


-1 


dt 




' a 




-1 


dt 


M 


M 


M 




M 




y'u-yl. 


' a 


-1 




^-4 






dt 



The first four columns of the above matrix Ha,5 are the same as that shown above in 
equation (2) and the fifth column comprises the partial derivative of each 
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transmitters pseudo-range, measurement with respect to the time-stamp. This partial 
derivative for the i-th transmitter may be expressed as: 



dt 



dt 



ix%(f,)-x:) + 



dt 



^l(xlih)-x:f +{yM-y:y +(z'M-z:f 



[1062] Similarly, the geometry matrices He^s for five measurements can be 
expressed as: 



SI 



yu-yst 2»-zsi 

~ yS2 ~ ^52 

M M 



M 

<-JCs„ y^-y^n z:-zL 
|ir'-5:| 



-1 



-1 



dt 

Ml 
dt 

M 



-I Ml 
dt 



Again, the first four columns of the above matrix H^^ are the same as that shown above 

in equation (3) and the fifth column comprises the partial derivative of each 

of 

transmitter's pseudo-range measurement with respect to the time-stamp. This partial 
derivative for the i-th transmitter may be expressed as: 



dxi 



dt 



dt 



dy's> 



t=tk 



dt 



iy%it,)-yl)+ 



dz% 



dt 



izM-zl) 



»=», 



4(^lih)-x:Y +{y'M-ylf +(z'M-zlf 



[1063] The partial derivatives (e.g., — ^, and of each satellite's 

at dt dt 

coordinate with respect to the measurement time can be derived from the Almanac or 
Ephemeris accordingly, since they described the satellite's trajectory as a function of 
time. For the above equations, is the estimate of the measurement time at the 
iteration. From the above equations, at each iteration, the satellite positions are 
recomputed based on the new value of the measurement time r, which may be obtained 
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from the measurement time eiror that is one of the unknowns whose value is reassessed 
at each iteration. 

[1064] The position determining calculations for the analog mode is described in 
further detail in U.S Patent Application Serial No. 09/280,337, entitled " METHOD 
AND APPARATUS FOR LOCATING GPS EQUIPPED WIRELESS DEVICES 
OPERATING IN ANALOG MODE," filed March 29, 1999, assigned to die assignee of 
the present application. 

[1065] In FIG. 2, processing circuits 212 and 222 may derive the initial coarse 
position estimate and the more accurate position estimate for the receiver device 110, 
respectively, by executing program instmctions stored within memories 214 and 224, 
respectively. Processing circuits 212 and 214 may each be implemented as a digital 
signal processor (DSP), an application specific integrated circuit (ASIC), a processor, a 
microprocessor, a controller, a microcontroller, a field programmable gate array 
(FPGA), a programmable logic device, other electronic unit, or any combination thereof 
designed to perform the functions described herein. Processing circuits 212 and 222 
may each further include memory 212 or 224 used to store program instmctions and 
data. 

[1066] The disclosed method and apparatus is provided to enable any person skilled 
in the art to make or use the present invention. Various modifications to these 
embodiments will be readily apparent to those skilled in the art, and the generic 
principles defined herein may be applied to other embodiments without departing from 
the spirit or scope of the invention. Thus, the present invention is not intended to be 
limited to the embodiments shown herein, but is to be accorded the broadest scope that 
is consistent with the disclosed principles and novel features. 



[10«7] WHATISCLAIMEBIS: 
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CLADVfS 



1. A method for determining the location of a device, comprising: 

2 (a) receiving an estimate of the location of the device, the estimate having been 

calculated using initial estimates of the position of a plurality of 
4 transmitters; 

(b) calculating a set of range estinoates using a revised estimate of the position of 
6 the device and revised estimates of the positions of the transmitters; 

(c) deriving an intermediate position estimate for the device based on the set of 
8 measurements and the initial position estimates for the transmitters; 

(d) computing an update vector based in part on the intermediate and revised 
1 0 position estimates for the device; and 

(e) updating the revised position estimate for the device based on the update 
12 vector. 

2. The method of claim 1, further comprising: 

2 repeating the calculating, deriving, computing, and updating a plurality of times. 

3. The method of claim 2, wherein the calculating, deriving, computing, and 
2 updating are repeated until the magnitude of the update vector is within a 

particular threshold. 

4. The method of claim 1, further comprising: 

2 (a) initializing the revised position estimate for the device to the coarse position 

estimate. 

5. The method of claim 1, wherein the transmitters are satellites in a Global 
2 Positioning Satellite (GPS) system. 

6. The method of claim 5, wherein the initial position estimates for the satellites are 
2 based on an Almanac. 



7. 

2 



The method of claim 5, wherein the revised position estimates for the satellites 
are based on an Ephemeris. 



wo 02/071095 



PCT/US02/06213 



20 



8. The method of claim 1, wheiein the revised position estimates for the 
2 transmitters are determined based on the time at which signals from the 

transmitters, which were used to derive the coarse position estimate for the 
4 device, were measured at the device. 

9. The method of claim 1, wherein the intermediate position estimate for the device 
2 is derived based on a least mean square O^JMS) computation. 

10. The method of claim 1, further comprising: 

2 (a) determining a first geometry matrix based on the initial position estimates for 

the transmitters and the intermediate position estimate for the device; and 

4 (b) determining a second geometry matrix based on the revised position 

estimates for the transmitters and the revised position estimate for the 

6 device, and 

(c) wherein the update vector is further based on the first and second geometry 

8 matrices. 

11. The method of claim 1, wherein the coarse and revised position estimates for the 
2 device are derived based on four measurements for four transmitters. 

12. The method of claim 11, wherein the measurements correspond to pseudo- 
2 ranges indicative of distances between the device and the transmitters. 

13. The method of claim 1, wherein the coarse and revised position estimates for the 
2 device are derived based on five or more measurements for five or more 

transmitters. 

14. The method of claim 13, wherein derivation of position estimate for the device is 
2 based on the five or more measurements to account for uncertainty in the time at 

which the measurements were made at the device. 



15. 



A method for determining the location of a device, comprising: 
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2 (a) receiving a coarse position estimate for the device derived based in part on 

initial estimates of the position of a plurality of transmitters; 
4 (b) receiving revised position estimates for the transmitters; 

(c) initializing a revised position estimate for the device; 

6 (d) coicputing an update vector based on the initial and revised position 

estimates for the device and the initial and revised position estimates for the 
8 transmitters; 

(e) updating the revised position estimate for the device based on the update 
10 vector, and 

(f ) repeating the computing and updating a plurality of times. 

16. A computer program product for determining the location of a device, 
2 comprising: 

(a) code for receiving a coarse position estimate for the device derived based in 
4 part on initial estimates of the position of a plurality of transmittals; 

(b) code for receiving revised position estimates for the transmitters; 
6 (c) code for initializing a revised position estimate for the device; 

(d) code for computing an update vector based on the initial and revised position 
8 estimates for the device and the initial and revised position estimates for the 

transmitters; 

10 (e) code for updating the revised position estimate for the device based on the 

update vector; 

12 (f ) code for repeating the computing and updating a plurality of times, and 

(g) a data storage medium configured to store the codes. 

17. A position identifying device comprising: 

2 (a) a communication port configured to receive a coarse position estimate for a 

receiver device derived based in part on initial estimates of the position of a 
4 plurality of transmitters; and 

(b) a processor operatively coupled to the communication port and configured 
6 to: 

(i) receive information indicative of revised position estimates for the 
8 transmitters, 

(ii) initialize a revised position estimate for the receiver device. 
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10 (iii) compute an update vector based on the initial and revised position 

estimates for the receiver device and the initial and revised position 
1 2 estimates for the transmitters, 

(iv) update the revised position estimate for the receiver cfevice based on 
14 the update vector, and 

(v) repeat the compute and update a plurality of times. 

18. The position identifying device of claim 17, wherein to compute the update 
2 vector, the processor is further configured to: 

(a) calculate a set of measurements based on the revised position estimate for the 
4 receiver device and the revised position estimates for the transmitters, 

(b) derive an intermediate position estimate for the receiver device based on the 
6 set of measurements and the initial position estimates for the transmitters, and 

(c) calculate the update vector based at least in part on the intermediate and 
8 revised position estimates for the receiver device. 

19. The position identifying device of claim 17, further comprising: 

2 (a) a memory operatively coupled to the processor and configured to store a 

plurality of sets of information indicative of position estimates for the 

4 transmitters at a plurality of time instances, and wherein the revised position 

estimates for the transmitters are derived from one of the sets of 

6 information. 

20. The position identifying device of claim 17, wherein the communication port is 
2 further configured to receive information indicative of a particular time at which 

signals from the transmitters were measured and used to derive the coarse 
4 position estimate for the receiver device. 
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